Interband proximity effect and nodes of superconducting gap in Sr2Ru04 
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The power-law temperature dependences of the specific heat, the nuclear relaxation rate, and 
the thermal conductivity suggest the presence of line nodes in the superconducting gap of S^RuCU. 
These recent experimental observations contradict the scenario of a nodeless (k x + ik y )-type super- 
conducting order parameter. We propose that interaction of superconducting order parameters on 
different sheets of the Fermi surface is a key to understanding the above discrepancy. A full gap 
exists in the active band, which drives the superconducting instability, while line nodes develop in 
passive bands by interband proximity effect. 

PACS numbers: 74.70.Pq, 74.20.Rp, 74.25.Bt 



The layered perovskite Sr2Ru04 with T c ~ 1.5 K 
[EJ is an example of an unconventional superconductor 
with non-s-wave Cooper pairing p|. The theoretical 
proposal of a spin-triplet p-wave order parameter 

A a| a(k) = (ia v a z ) a f3 d(k), d(k) cx (k x + ik y ) is supported 
by experimental observations of a temperature indepen- 
dent Knight shift for H _L c and an increased muon 
spin-relaxation below T c [Q. Such an axial gap function 
has a nonvanishing amplitude on the cylindrical quasi- 
two-dimensional Fermi surface of Sr 2 Ru04 [^-|iot|. This 
property favors the axial state as a natural choice in a 
weak-coupling theory, which generally supports nodeless 
solutions H . Recent experimental data collected on high 
quality samples, however, seem to invalidate the above 
conclusion. The power-law temperature dependences as 
T -> found for the specific heat, C(T) cx T 2 
the NQR relaxation rate, Tf 1 cx T 3 the thermal 

conductivity, k(T) cx T 2 |14}l5|], the penetration depth 
Jl6| , and the ultrasonic attenuation |l7[] point to lines of 
zeros in the superconducting gap and, thus, question the 
consistency of the whole picture. 

There have been several theoretical attempts to resolve 
this controversy |lg|-|2l]] . Most suggest replacing the axial 
p-wave order parameter d p (k) cx (k x + ik y ) by a suitable 
/-wave gap: rf/(k) cx {k x + ik y )g(k), where the even par- 
ity function g (k) is chosen to have zeros, e.g. k x k y or 
(k 2 — k 2 ) |19|-pi|. There is no clear microscopic mecha- 
nism for such an /-wave instability. More importantly, 
nodes in an /-wave gap are only marginally stable, i.e. 
they disappear if all symmetry allowed harmonics are in- 
cluded in the expansion of the gap function. For example, 
for g(k) = k x k y one finds: 



d(k) = T]i (k x + ik y )k x ky + ir\i(k x - ik y ) 



(1) 



where r/i and 772 are real. Both terms in Eq. ([!]) trans- 
form in the same way under operations of the symmetry 
group of the superconducting state. In particular, both 
harmonics are symmetric with respect to a four-fold rota- 
tion e l7r / 2 C4 and time-reversal in combination with a re- 



flection in the (100) plane, e lv Ta x . Therefore, the /- and 
the p-wave harmonics are symmetry indistinguishable in 
tetragonal crystals and mix with each other producing a 
finite gap |A| m ; n ~ 7/2- Very recently, Izawa et al. have 
measured the basal plane anisotropy of the thermal con- 
ductivity k(6) in Sr2RuC>4 at finite magnetic fields [flB] . 
Their results also discard an /-wave gap together with a 
so-called 'anisotropic' p-wave state |18|] as possible can- 
didates to explain line of nodes in S^RuC^: all these 
superconducting states have a substantial anisotropy of 
k(9) in magnetic field determined by in-plane node struc- 
ture, whereas experimentally the basal plane anisotropy 
is much smaller. Izawa et al. |1| suggest instead hori- 
zontal lines of nodes in the superconducting gap. 

An opposite conclusion has been reached by Lupien et 
al. |l7]] from the anisotropy of the ultrasonic absorption. 
However, the measured anisotropy appears mainly in the 
absolute magnitude of the attenuation and not in the ex- 
ponent of the temperature power law. For this reason it 
is not clear that these results are in conflict with hori- 
zontal lines of zeros and, as Lupien et al. stress, detailed 
calculations based on the actual electronic structure are 
necessary for a definitive interpretation of their results. 

Here, we propose a mechanism for the formation of 
horizontal line nodes in the superconducting gap of 
Sr2RuC>4. The Fermi energy crosses three bands, deter- 
mined by the d xy (7-sheet of the Fermi surface) and the 
hybridized d xz and d yz (a- and /3-sheets) orbitals of Ru 
|§-|Io|. Magnetic fluctuations, responsible for anisotropic 
Cooper pairing H|2^1 , have significant orbital dependence 
p3|. Therefore, the intrinsic temperature of the super- 
conducting instability should vary from band to band 
with one sheet being the active source for superconduct- 
ing instability and the others being the passive sheets. In 
reality, interband scattering of Cooper pairs, or proxim- 
ity effect in the momentum space, will induce the super- 
conducting gap simultaneously on all parts of the Fermi 
surface. Such interband scattering is generally a strong 
effect, which allows one to treat numerous multiband su- 
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perconductors by an effective single band Fermi surface. 
Agterberg and co-workers B have argued that Sr2Ru04 
is different: a direct in-plane scattering of the p-wave 
Cooper pairs between bands is significantly suppressed 
by the orbital symmetry. Therefore they conclude, one 
or two bands develop only tiny superconducting gaps, 
which show up at intermediate temperatures as a residual 
density of states with a subsequent crossover as T — > 
to a full gap behavior. In this Letter we study addi- 
tional interlayer contributions to interband scattering of 
Cooper pairs, which become important when direct in- 
plane scattering is suppressed. We find that a nodeless 
axial order parameter d p (k) oc (k x + ik y ) in the active 
band can induce superconducting gaps with zeros in the 
passive bands: d' p (k) oc (k x + ik y ) cos(fc z /2). Thus, circu- 
lar nodes of the superconducting gap develop about the 
c axis on one or two of the three Fermi surface sheets. 
This model of weakly-coupled superconducting order pa- 
rameters in different bands fits well C(T) in zero field 
and helps to explain the observed field behavior of 
the specific heat [ fijl . 

We start with a general two-particle interaction: 

V = J dTdr'U(T,T')i(}\r)^(r')tp(r')tp(r) , (2) 

where for simplicity we omit all spin indices assuming a 
fixed spin structure of the triplet order parameter. The 
band representation for the effective interaction V is ob- 
tained by (i) expanding the field operators ip( r ) m terms 
of the band operators: V( r ) = J2i k Wk(r)c;k (I is a band 
index) and (ii) representing the Bloch function <piu{y) in 
a given band as a lattice sum over the Wannier function 
of the Ru orbitals: <p lu (r) = £„ e ' kR "<M r ~ R «)- The 
interaction in the Cooper channel is 

V = \ E V lv {k,k!)c\^c v ^c V u , (3) 

ll',kk' 

where the scattering vertex is given by 
Vi V (k,k')=[drdr'U{r,r') e - ik ( R "- R ™' V k '( R ™- R -') 

nn'mm' 

x#(r - R„)<^(r' - R„-)<Mr' - Rm')#(r - Rm). (4) 

Following Ref. 0] we now consider interband scattering 
processes (/ ^ in the tight-binding approximation, i.e. 
assume that the Wannier functions are well localized and, 
therefore, the main contribution to Vw (k, k') comes from 
a few neighboring sites. The largest on-site contribution 
(Rn = R n ' = R m = R m ') is independent of k and k'. It 
causes coupling only between conventional s-wave order 
parameters in two bands. The coupling of the p-wave 
order parameters appears first in the sum (^) for R m = 
R„, R m ' =R„', and (R„< - R„) = Si = ±ax(y) (a is 
lattice constant): 



Vj p (k,k')= sinfciasin/c^a 

i— x 1 y 

x Jdrdr f U(T,i>)4>t(T)MT)ti(''-*i)M''-ti)- ( 5 ) 

This direct in-plane scattering of the Cooper pairs in- 
duces the same nodeless superconducting gap 

di(k) oc (sin k x a + i sin k y a) (6) 

on all sheets of the Fermi surface. If we now approx- 
imate in a tight-binding spirit U(r,r') « U(di), then 
the double integral in Eq. (g) factorizes in a product of 
two spatial integrals. Each of these integrals vanishes 
separately for I = 7 and l' = a,(3 because the orbitals 
from these bands have different parity with respect to 
a z . Thus, it is essential to keep spatial dependence of 
U(r, r'). For a Coulomb-type interaction the off-diagonal 
matrix element (/ ^ I') in Eq. (||) has a dipolar reduc- 
tion {b/a) 2 ~ 0.02 (J] compared to the diagonal matrix 
elements (I — I'), where & is a characteristic spatial ex- 
tension of the Wannier functions. In reality, the matrix 
element will be reduced even further because d xz and d yz 
orbitals mix in a- and /3-bands only in the close vicinity of 
the Brillouin zone diagonals. Away from these directions 
there is an extra approximate symmetry (T x (y)i which in- 
troduces an effective quadrupolar reduction of the matrix 
element in Eq. (||). Thus, the direct interaction of the p- 
wave order parameters between 7- and a- or /3-bands is 
significantly reduced and the amplitude of the type-I gap 
(g) in passive bands is much smaller than in the active 
band ||. 

The next contribution to the interband scattering of 
the p-wave pairs in Eq. (^J) comes from interlayer terms 
with R m = R„ and R n / = R„ ± ax(y ) , R m / = R„ ± |x ± 
|y ± |z on a bet lattice of Ru-atoms. Summing over all 
contributing sites the p-wave gap of the type (||) in the 
7-band induces 

a . k v a k x a\ k z c 

(folk) oc sin cos — - — I- i sin — — cos cos 

K ' V 2 2 2 2 / 2 

(7) 

in a- and /3-bands and vice versa. Existence of the type- 
II p-wave gap, but in all bands simultaneously, has been 
conjectured by Hasegawa et al. [^9|. They, however, 
based their suggestion on an (unjustified) assumption of 
a repulsive interaction between electrons in a single Ru- 
O plane and an attraction only for electrons in adjacent 
layers. The type-II superconducting gap d 2 (k) has circu- 
lar line nodes at k z — ±7r/c. Importantly, they are stable 
with respect to an admixture of a small amount of the 
type-I gap, which only shifts the position of zeros along 
k z . The two gaps are mixed with a real phase, as re- 
quired by the time-reversal symmetry (Ta x ), and nodes 
of di + d 2 disappear only for |di| roax > |d 2 |max- 
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A reliable estimate for the strength of the scatter- 
ing vertices corresponding to the two types of induced 
gaps @ and (fjj) must be done in a truly microscopic 
treatment of the Fermi liquid state in S^RuCj, which 
is beyond the scope of our work. Note that scatter- 
ing processes contributing to Eq. (Q) have one sym- 
metry cancellation factor less than Eq. (||), but in- 
stead they are reduced by the small overlap of the 
orbitals in adjacent layers. Information on interlayer 
overlap can be obtained by analyzing the results of 
the high-precision de Haas-van Alphen measurements 
JlO| |, which determined corrugation of the Fermi surface 
cylinders along k z . Bergemann et al, flo[| found the 
strongest corrugation in the /3-sheet of the Fermi sur- 
face: Akp^f} ~ cos(k z c/2) and a much weaker corruga- 
tion of the 7-sheet: Akp n ~ cos(fczc). This different 
periodicity naturally appears in a tight-binding model, 
since d xz (d yz ) orbitals have a direct interlayer overlap 
t'j_, leading to a diagonal contribution to the band en- 
ergy: 8t' ± cos(k x a/2) cos(k y a/2) cos(k z c/2). On the other 
hand, planar d xy orbitals do not hybridize across the lay- 
ers. However, they can hybridize with d xz (d yz ) orbitals 
in adjacent planes with the matrix element t" , which con- 
tributes 8t'^cos(k x (y)a/2)sm(k y ( x )a/2)sm(k z c/2) to the 
off-diagonal kinetic energy. As a result, the corrugation 
of the /3-cylinder is a first-order effect in t' ± , whereas the 
corrugation of the 7-sheet has a much weaker second- 
order contribution ~ cos(k z c/2)] 2 /t. Comparing the 
experimental corrugations we find: t'j_ ~ — 1 meV and 
somewhat larger magnitude for <" ~ 3 meV. Since the 
interlayer hopping amplitude is proportional to the or- 
bital overlap in Eq. (^) , we argue that it is quite possible 
to have a comparable amplitude for the two types of the 
p-wave gaps induced by interband proximity effect in pas- 
sive bands. 

We consider now the effect of line nodes in passive 
bands on thermodynamic properties of the superconduct- 
ing state. In particular, we calculate the specific heat 
C(T). Since the Cooper pair scattering between the a- 
and the /3-sheets is not small, we adopt an effective two 
band model for Sr2Ru04 and split the total density of 
states at the Fermi level according to AT 01 : N 02 — N 0l : 
(No a + Nop) — 0.57 : 0.43 based on the de Haas-van 
Alphen measurements ||. We also assume that the ac- 
tive sheet for the superconducting instability is the 7- 
sheet. Our main motivation for this assumption comes 
from comparison with the experimental data below. We 
adopt a weak-coupling approach and parameterize the 
pairing potential in the two bands by three parameters: 

y n (k,k') = - 5l /(k)/(kQ, y 22 (k,k') = - 52 /(k)/(k'), 

and Vi2(k, k') = — g 3 /(k)/(k'), where we choose for sim- 
plicity /(k) = k/fcp and /(k) = V2(k/k F ) cos(fc z /2), i.e. 
we presume that only interlayer processes contribute to 
the interband scattering of the Cooper pairs. The pair- 
ing interaction in the active band is attractive (g\ > 0), 
while interaction constants in the passive band ((72) and 



between the bands (33) can have arbitrary sign. Solv- 
ing the system of the two gap equations numerically we 
determine the specific heat from 



C(T) = 2Y,E lk - 



df(E, 



Ik) 



dT 



(8) 



where -E/k is the quasiparticle excitation energy 
\Af k + A^(k) (we consider only unitary triplet states) 
and /(-Eik) is the corresponding Fermi distribution. 
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FIG. 1. Temperature dependence of the normalized spe- 
cific heat. The upper panel: the two-band model results for 
various choices of the interaction parameters. Curves #1-3 
correspond to (72/ffi = 0.85 and = 0.01,0.07,0.2, re- 

spectively. The curve #4 is for (72/(71 = 0.1 and (73/(71 = 0.07. 
The lower panel: circles are experimental data for S^RuOi 
[11]. One-band results are shown for anisotropic gap with line 
nodes (dashed line) and for isotropic gap (dot-dashed line). 
Solid line is the two-band model fit with (72 — 0.85gi and 
S3 =0.07 ff i. 

We present in the upper panel of Fig. 1 the specific 
heat for different choices of the coupling constants in the 
two-band model. The calculations have been done for a 
'typical' weak-coupling magnitude of gi = 0.4 (in units 
of the inverse total density of states) varying the two 
other parameters. The first three curves correspond to 
weak (73 = O.Olgi (#1), intermediate g 3 = 0.07gi (#2), 
and strong g 3 — 0.2gi (#3) interband scattering of the 
Cooper pairs, keeping in all cases (72 = 0.85<?i. Such a 
moderate change in the coupling constants between the 
two bands A2/A1 = (72-^02/51-^01 = 0.64 results in an or- 
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der of magnitude difference in their bare transition tem- 
peratures: T c ijT c \ — 0.086 in the weak-coupling theory. 
For the weaker interband coupling the heat capacity de- 
velops a second peak, which reflects a nonzero bare tran- 
sition temperature in the passive band. For the stronger 
interband coupling the two gaps are tightly bound to 
each other and we return to an effective one-band behav- 
ior. The curve #4 is an example of a shoulder in the 
temperature dependence of C/T, which arises for a re- 
duced pairing interaction in the passive band 52 = O.I51 
when we keep the same scattering vertex 53 = 0.07gi as 
for the curve #2. 

In the lower panel of Fig. 1 we present our fit to the ex- 
perimental data of NishiZaki et al. jL2|, which coincides 
with the curve #2 above, and also show the results of one- 
band models with constant and anisotropic gaps. Under 
condition of a reduced interband coupling the heat capac- 
ity jump at T c is entirely determined by the active band. 
Remarkable agreement between the experimental mag- 
nitude of the jump and AC/C n ~ 1A3N 7 /N Q = 0.82, 
which is obtained from the standard BCS result by its 
renormalization on the partial density of states of the 7- 
sheet, is a direct confirmation of our choice of the active 
band. Though, it is impossible to fix all three parameters 
of the two-band model uniquely, this model can naturally 
explain a clear convex shape of the experimental data for 
C /T at low temperatures by choosing an intermediate 
strength of the interband scattering matrix element. In 
contrast, a one-band model with anisotropic gap or the 
two-band model with a strong interband scattering pre- 
dict a concave shape for C/T. From our fit we cannot also 
exclude a possibility of a small but finite |A m ; n | in the 
passive bands, which appears if in-plane scattering am- 
plitude slightly exceeds the interplane contribution (|?]). 

The field dependence of the residual density of states 
at low temperatures suggests another argument in favor 
of the multiband scenario. Small magnetic fields (<C H C "i) 
quickly restore about 40% of the total density of states 
In our model this corresponds to the behavior ex- 
hibited by the curve #4, the upper panel in Fig. 1, in 
the case of temperature effects. Such a new feature in 
external field arises because of an additional suppression 
of superconductivity in the passive a- and /3-bands for 
H _L c. Stronger c-axis dispersion in these bands leads 
to a larger coherence length £ c and an extra reduction 
of the bare H^- This effect should disappear for H || c 
because of similar values of the in-plane Fermi velocities, 
which also agrees with the experiment jl2]. It would be 
also interesting to reinvestigate the impurity effect on 
the residual density of states. Such an analysis has been 
done previously for the two-band model with a constant 
gap amplitude in the passive band pq| . Line nodes can 
modify the expected behavior and produce a gapless su- 
perconducting state in the passive bands. 

In summary, we have shown that circular horizontal 
line nodes in the superconducting gap of Sr2Ru04 ap- 



pear due to weak and anisotropic interband proximity 
effect. This effect is a consequence of (i) non s-wave sym- 
metry of the Cooper pairs [conventional superconductors 
have generally a strong isotropic interband coupling dom- 
inated by on-site term in Eq. (|j)] and (ii) specific sym- 
metry of the Ru orbitals, which give extra suppression 
of the matrix element in Eq. (|^). Further experimen- 
tal tests of our scenario should include studying effects 
of pressure, which can modify the strength of interlayer 
scattering amplitude for the Cooper pairs. 
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supported by Swiss National Fund. 
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